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1.0  Overview 

1.1  Summary  of  Research  -  Statement  of  Work 

(1)  Develop  generalized  field  transforms  that  provide  the  basis 

for  complete  expansions  for  the  electromagnetic  fields  in  rough  surfaces. 
The  generalized  transforms  will  be  used  to  reduce  Maxwell's  electro¬ 
magnetic  equations  in  conjunction  with  exact  boundary  conditions  into 
sets  of  first  order  coupled  differential  equations  for  forward  and 
backward  travelling  wave  amplitudes. 

The  elements  of  the  full  wave  method,  its  principal  properties  and 
its  relationships  to  the  physical  optics  and  perturbation  solutions 
are  presented  in  the  RADC-TR-82-336  Annual  Report  Section  1.1  (1  March 
1981-28  February  1982)  The  generaiiZecj  field  transforms  and 

the  coupled  differential  equations  for  the  forward  and  backward 
travelling  wave  amplitudes  are  presented  in  the  Second  Annual  Report 
Section  1.1  (1  March  1982-28  February  1983). (26) 

(2)  The  iterative  solutions  for  the  vertically  and  horizontally 
polarized  scattered  radiation  fields  when  the  slope  of  the  rough 
surface  is  small  will  be  obtained  through  the  use  of  the  steepest 
descent  method. 

The  second  order  iterative  solutions  for  the  scattered  radiation 
fields  are  presented  in  the  Second  Annual  Report 

(3)  Utilize  the  results  of  (2)  extend  the  studies  to  obtain  full 
wave  solutions  for  scattered  radiation  fields  from  two  dimensionally 
rough  surfaces  with  arbitrary  slopes. 

The  full  wave  solutions  for  the  scattered  radiation  fields  from 
two  dimensionally  rough  surfaces  are  presented  in  the  Second  Annual 
Report,  Section  (1.1) v  .  This  work  has  also  been  applied  to  problems 

The  publications  cited  in  this  section  (1)  through  (27)  are  listed 
in  Section  1.2  of  this  report. 
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of  depolarization  of  the  scattered  radiation  field  by  conducting 

(3) 

objects  of  irregular  shape  above  rough  land  and  sea  .  The  full 
wave  approach  was  also  used  to  determine  the  scattering  and 

depolarization  of  radio  waves  in  irregular  spheroidal  struc- 

(5), (6), (12). 
tures  *  * 

(4)  The  full  wave  solution  obtained  in  (3)  will  be  compared 
with  the  perturbational  and  physical  optics  solutions. 

The  comparison  of  the  full  wave  solutions  with  the  perturba¬ 
tion  and  physical  optics  solutions  is  presented  in  the  Interim 
Technical  Report  RADC-TR-81-203  (August  1981) .  This  work  has  also 
been  used  to  resolve  the  discrepancies  between  different  physical 

/«\  /  ry  y  \ 

optics  solutions  for  rough  surface  scattering  ' ’ '  % 

(5)  The  normalized  backscatter  cross  section  for  an  arbitrary 
rough  surface  will  be  obtained  from  the  results  of  (2),  (3),  and  (4). 

The  normalized  like  and  cross  polarized  backscatter  cross 
sections  are  derived  for  arbitrary  rough  surfaces  ^  ^  . 

(6)  The  polarization  dependence  of  the  backscatter  cross 
section  grazing  incidence  and  incidence  at  pseudo-Brewster 
angle  with  respect  to  reference  plane  will  be  studied. 

The  polarization  dependence  of  the  scattering  cross  sections 

are  derived  for  all  angles  of  incidence^\  These  solutions 

(9) 

are  compared  with  the  perturbed-physical  optics  technique 

(7)  The  study  will  emphasize  the  relationship  of  scattered  fields 
to  the  statistics  of  the  height  and  slope  of  the  rough  surface  which 
are  given  by  spectral  density  of  the  surface  height,  the  correlation 
function  and  the  correlation  distance.  The  study  will  take  into 
account  the  absorption  of  electromagnetic  energy  by  the  medium. 
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The  scattering  cross  sections  have  been  computed  for  rough 

surfaces  that  are  specified  through  their  surface  height  spectral 

(8) 

density  functions  of  the  surface  height  autocorrelation  function. 

This  work  is  based  on  a  two-scale  model  of  the  composite  rough 

surface.  Thus  the  problem  of  specifying  the  wavenumber  where 

spectral  splitting  occurs  is  considered  in  detail. 

A  uniform  full  wave  approach  that  does  not  adopt  hybrid  two- 

(13) 

scale  models  of  rough  surfaces  has  also  been  developed  recently 

Details  of  this  formulation  are  presented  in  section  2  of  this 
(27) 

report  .  In  this  work,  the  surface  is  characterized  by  a  complex 
dielectric  coefficient  to  account  for  absorption  of  electromagnetic 
energy  by  the  medium. 

Scattering  by  Non-Gaussian  rough  surfaces  has  also  been 

presented  in  the  RADC-TR-83-47  Interim  Technical  report  (March  1983) . 

In  this  work  the  shadow  function  and  the  backscatter  cross  sections 

have  been  computed  for  Non-Gaussian  rough  surfaces  for  which 

decorrelation  implies  statistical  independence  . 

In  addition  to  publications  in  technical  and  scientific  journals 

through  and  the  Interim,  Annual  and  Final  Reports 

(27) 

through  ,  the  principal  investigator  has  presented  invited  and 

(14) 

contributed  papers  at  several  international  conferences  through 

(22) 

.  Ten  quarterly  reports  were  also  submitted  during  this 
period . 
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1.2  List  of  Publications  by  the  Principal  Investigator  During  the 
Period  of  the  Contract 

1.2.1  Publications  in  Technical  Journals 

1.  "Scattering  Cross  Sections  for  Composite  Surfaces  with  Large 
Mean  Square  Slopes — Full  Wave  Analysis,"  International  Journal 
of  Remote  Sensing,  Vol.  3,  No.  3,  pp.  327-337,  1982. 

2.  "Resolution  of  the  Discrepancies  Between  Different  Physical 
Optics  Solutions  for  Rough  Surface  Scattering,"  Journal  of  the 
Acoustical  Society  of  America.  Vol.  71,  No.  6,  June  1982, 

pp.  1352-1358. 

3.  "Depolarization  of  the  Scattered  Radiation  Field  by  Conducting 
Objects  of  Irregular  Shape  Above  Rough  Land  and  Sea — Full  Wave 
Solutions,"  Radio  Science,  Vol.  17,  No.  5,  September-October 
1982,  pp.  1055-1066. 

4.  "Scattering  and  Depolarization  by  Rough  Surfaces:  Full  Wave 
Approach,"  Proceedings  of  the  SPIE  International  Technical 
Symposium  of  the  International  Society  of  Optical  Engineering, 
Vol.  358,  No.  28,  pp.  1-14,  August  1982. 

5.  "Propagation  of  Vertically  and  Horizontally  Polarized  Waves 
Excited  by  Distributions  of  Electric  and  Magnetic  Sources  in 
Irregular  Stratified  Spheroidal  Structures  of  Finite 
Conductivity  Generalized  Field  Transforms.  Canadian  Journal 
of  Physics.  Vo.  61,  No.  1,  pp.  113-127,  1983. 

6.  "Scattering  and  Depolarization  of  Electromagnetic  Waves  in 
Irregular  Stratified  Spheroidal  Structures  of  Finite 
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Conductivity — Full  Wave  Analysis.  Canadian  Journal  of  Physics. 

Vol.  61,  No.  1,  pp.  128-139,  1983. 

"Scattering  Cross  Sections  for  Composite  Surfaces  that  Cannot 
be  Treated  as  Perturbed — Physical  Optics  Problems."  Radio 
Science .  Vol.  18,  No.  2,  pp.  129-137,  April  1983. 

"Computations  of  Scattering  Cross  Sections  for  Composite 
Surfaces  and  the  Specification  of  the  Wavenumber  Where  Spectral 
Splitting  Occurs,"  IEEE  Transactions  on  Antennas  and 
Propagation — in  press. 

"Comparison  of  Backscatter  Cross  Sections  for  Composite  Rough 
Surfaces  with  Different  Mean  Square  Slopes,"  International 
Journal  of  Remote  Sensing — in  press. 

"Shadowing  by  Non-Gaussian  Rough  Surfaces  for  Which 
Decorrelation  Implies  Statistical  Independence."  Radio  Science, 

Vol.  18,  No.  4,  pp.  566-572,  July-August  1983. 

"Scattering  Cross  Sections  for  Composite  Models  of  Non-Gaussian 
Rough  Surfaces  for  Which  Decorrelation  Implies  Statistical 
Independence,"  IEEE  Transactions  on  Antennas  and  Propagation — in  press. 
"Excitation,  Scattering  and  Depolarization  of  Electromagnetic 
Waves  in  Irregular  Stratified  Structures  of  Finite  Conductivity 
— Full  Wave  Analysis,"  Proceedings  of  the  International  Union 
of  Radio  Science  URSI  Conference  at  Ciudad  Universitaria, 

Madrid,  August  1983,  in  press. 

"Scattering  Cross  Sections  for  Composite  Rough  Surfaces  Using 
the  Unified  Full  Wave  Approach."  Submitted  for  review. 


•*«  ^ • 


1.2.2  Abstracts  and  Summaries  of  Papers  Presented  at  International 


Conferences 

14.  International  IEEE/APS  Symposium  and  National  Radio  Science 
Meeting,  Los  Angeles,  June  15-19,  1981,  "Scattering  Cross 
Sections  for  Composite  Random  Surfaces — Full  Wave  Solutions." 

15.  International  Union  of  Radio  Science  Symposium  on  Mathematical 
Models  in  Radio  Propagation  in  Conjunction  with  the  XXI 
General  Assembly  of  URSI.  Invited  Review  Paper  on  "Depolari¬ 
zation  of  the  Scattered  Radiation  Fields  by  Conducting  Objects 
of  Irregular  Shape  Above  Rough  Land  and  Sea — Full  Wave 
Solutions,"  Washington.  D.  C. .  August  10-19,  1981, 

16.  International  Union  of  Radio  Science  Open  Symposium  on  Remote 
Sensing  August  10-19,  1981,  Washington,  D.  C.,  "Scattering 
Cross  Sections  for  Random  Two  Dimensionally  Rough  Surfaces; 
Perturbation  Physical  Optics  and  Full  Wave  Solutions." 

17.  International  Union  of  Radio  Science  (URSI)  Symposium  at  the 
University  of  Colorado,  Jan.  13-15,  1982,  "Resolution  of  the 
Discrepancies  Between  Different  Physical  Optics  Solutions  for 
Rough  Surface  Scattering." 

18.  Joint  International  IEEE/APS  and  National  Radio  Science 
Meeting  at  the  University  of  New  Mexico,  May  24-28,  1982, 

"Scattering  Cross  Sections  for  Composite  Surfaces  with  Large 
Mean  Square  Slopes — Full  Wave  Analysis." 

19.  The  SPIE  26th  Annual  International  Technical  Symposium  of  the 
International  Society  for  Optical  Engineering,  San  Diego,  Calif.. 
Aug.  23-27,  1982.  "Scattering  and  Depolarization  by  Rough  Surfaces. 


20. 
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International  IEEE/APS  Symposium  and  National  Radio  Science 
Meeting  at  the  University  of  Houston,  Texas,  May  23-26,  1983, 
"Rough  Surface  Scattering  that  Cannot  be  Analyzed  Perturbed- 
Physical  Optics  Approaches." 

21.  International  IEEE/APS  Symposium  and  National  Radio  Science 
Meeting  at  the  University  of  Houston,  Texas,  May  23-26,  1983, 
"Scattering  Cross  Sections  for  Composite  Surfaces  and  the 
Wavenumber  Where  Spectral  Splitting  Occurs." 

22.  International  Union  of  Radio  Science  (URSI)  Meeting  at  Ciudad 
Universitaria,  Madrid,  August  23-26,  1983,  "Excitation, 
Scattering  and  Depolarization  of  Electromagnetic  Waves  in 
Irregular  Stratified  Spheroidal  Structures  of  Finite 
Conductivity — Full  Wave  Analysis." 

1.2,3  Annual  Technical  Reports  and  Interim  Technical  Reports 

23.  "Full  Wave  Analysis  for  Scattering  Cross  Sections,  Part  1: 
Random  Rough  Surfaces  and  Part  2:  Composite  Surfaces," 

U.  S.  Air  Force  Interim  Report,  August  1981.  RADC-TR-81-204 . 

(60  pages) . 

24.  "Scattering  and  Depolarization  of  Electromagnetic  Waves-Full 
Wave  Solutions,"  U.  S.  Air  Force  Annual  Report,  February 
1982  (90  pages) . 

25.  "Scattering  Cross  Sections  for  Composite  Models  of  Non- 
Gaussian  Surfaces  for  Which  Decorrelation  Implies  Statistical 
Independence,"  U.  S.  Air  Force,  Interim  Report,  November 
1982  (45  pages). 


"Scattering  and  Depolarization  of  Electromagnetic  Waves-Full 
Wave  Solutions,"  U.  S.  Air  Force  Annual  Report,  March  1,  1982 
February  28,  1983  (130  pages). 

"Scattering  and  Depolarization  of  Electromagnetic  Waves-Full 
Waves  Solutions,"  U.  S.  Air  Force  Final  Report,  March  1,  1981 
August  31,  1983  (119  pages). 


2.0  SCATTERING  CROSS  SECTIONS  FOR  COMPOSITE  ROUGH  SURFACES 


USING  THE  UNIFIED  FULL  WAVE  APPROACH 
2.1  Background 

The  full  wave  approach  is  used  to  derive  a  unified  formulation 
for  the  like  and  cross  polarized  scattering  cross  section  of 
composite  rough  surfaces  for  all  angles  of  incidence.  Earlier 
solutions  for  electromagnetic  scattering  by  composite  random  rough 
surfaces  are  based  on  two -scale  models  of  the  rough  surface.  Thus 
on  applying  a  hybrid  approach  Physical  Optics  theory  is  used  to 
account  for  specular  scattering  associated  with  a  filtered  surface 
(consisting  of  the  large  scale  spectral  components  of  the  surface) 
while  perturbation  theory  is  used  to  account  for  Bragg  scattering 
associated  with  the  surface  consisting  of  the  small  scale  spectral 
components . 

Since  the  full  wave  approach  accounts  for  both  specular  point 
scattering  and  Bragg  scattering  in  a  self  consistent  manner  the  two 
scale  model  of  the  rough  surface  is  not  adopted  in  this  work. 

These  unified  full  wave  solutions  are  compared  with  the  earlier 
solutions  and  the  simplifying  assumptions  that  are  common  to  all 
the  earlier  solutions  are  examined.  It  is  shown  that  while  the 
full  wave  solutions  for  the  like  polarized  scattering  cross 
sections  based  on  the  two-scale  model  are  in  reasonably  good 
agreement  with  the  unified  full  wave  solutions,  the  two  solutions 
for  the  cross  polarized  cross  sections  differ  very  significantly. 
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2.2  Discussion 

Traditionally  physical  optics  and  perturbation  theories  have 

been  used  to  derive  the  like  and  cross  polarized  scattering  cross 

sections  for  composite  random  rough  surfaces  (Beckmann  1963,  Rice 

1951).  To  this  end  two-scale  models  have  been  adopted  and  the 

rough  surfaces  are  regarded  as  small  scale  surface  perturbations 

that  are  superimposed  on  large  scale,  filtered  surfaces  (Wright 

1968,  Valenzuela  1968,  Barrick  and  Peake  1968).  Thus  the  scattering 

cross  sections  are  expressed  as  sums  of  two  cross  sections.  The 

first  accounts  for  specular  point  scattering.  It  is  given  by  the 

physical  optics  cross  section  for  the  filtered  surface  consisting 

of  the  large  scale  spectral  components.  The  second  accounts  for 

Bragg  scattering.  It  is  given  by  the  perturbation  cross  section 

for  the  surface  consisting  of  the  small  scale  spectral  components 

that  ride  on  the  filtered  surface. 

On  applying  the  perturbed-physical  optics  approaches  it  is 

necessary  to  specify  the  wavenumber  k^  where  spectral  splitting 

is  assumed  to  occur  between  the  large  and  small  scale  spectral 

components  of  the  rough  surface.  Thus  Brown, (1978)  who  applied 

a  combination  of  Burrows'  perturbation  theory  (1967)  and  physical 

optics  (Beckmann  1968),  to  obtain  the  scattering  cross  sections 

for  perfectly  conducting  random  rough  surfaces,  specified  k,  on 

d 

the  basis  of  the  characteristics  of  the  small  scale  surface 
2  2 

(8  =  4k  <h  >  ■  0.1,  where  k  is  the  electromagnetic  wavenumber  and 
os  o 

2 

<h  >  in  the  mean  square  height  of  the  small  scale  surface) . 
s 


However,  using  the  approaches  of  Hagfors  (1966)  and  Tyler  (1976) 
the  specification  of  k^  is  assumed  to  be  based  on  the  characteristics 
of  the  large  scale  surface.  In  general  the  restrictions  on  both 
the  large  and  small  scale  surfaces  cannot  be  satisfied  simultaneously 
and  using  the  perturbed-physical  optics  approaches  the  evaluation 
of  the  scattering  cross  sections  critically  depend  on  the  specifi¬ 
cation  of  kd  (Brown  1978). 

More  recently  the  full  wave  approach  has  been  used  to  determine 
the  scattering  cross  sections  for  composite  random  rough  surfaces 
of  finite  conductivity  (Bahar  1981b,  Bahar  and  Barrick  1983) . 

Since  the  full  wave  solutions  account  for  Bragg  scattering  and 
specular  point  scattering  in  a  self-consistent  manner,  it  is  not 
necessary  to  decompose  the  surface  into  two  surfaces  with  small 
and  large  roughness  scales.  However,  when  such  a  decomposition 
is  feasible,  the  full  wave  solutions  for  the  scattering  cross  sections 
can  be  expressed  in  terms  of  a  weighted  sum  of  two  cross  sections 
(Bahar  1981b,  Bahar  and  Barrick  1983)  .  Thus  on  adopting  a  two-scale 
model,  the  full  wave  solution  resolves  the  discrepancies  between 
Valenzuela's  (1968)  solution  (mostly  based  on  physical  considerations) 
and  Brown's  solution  (1978).  Furthermore,  in  an  attempt  to  draw 
more  definite  conclusions  regarding  the  choice  of  k^,  it  was  varied 
over  a  wide  range  of  values  (Bahar  et  al  1983) .  It  was  shown  that 
while,  as  expected,  the  cross  sections  associated  with  the  individual 
large  and  small  scale  surfaces  critically  depend  upon  the  choice  of 
k^,  the  weighted  sum  of  the  like  polarized  cross  sections  remain 
practically  insensitive  to  variations  in  k^  for  8  >  1.0.  Thus, 
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provided  that  the  large  scale  surface  satisfies  the  radii  of 
curvature  criteria  (associated  with  the  Kirchhoff  approximations 
for  the  surface  fields)  and  the  condition  for  deep  phase  modulation, 
the  full  wave  solutions  for  the  like  polarized  scattering  cross 
sections  based  on  the  two- scale  model  are  practically  independent 
of  the  specified  value  of  k^.  However,  on  applying  the  full  wave 
approach  to  evaluate  the  like  and  cross  polarized  scattering  cross 
sections  for  two- scale  models  of  composite  rough  surfaces,  several 
assumptions  were  made  to  facilitate  the  computations.  The  first 
assumption  was  that  the  large  and  small  scale  surfaces  were  statistically 
independent  (Brown  1978) .  It  would  seem  reasonable  to  make  such 
an  assumption  if  the  two  surfaces  are  results  of  independent  processes. 
This  would  be  the  case,  for  example,  if  the  small  scale  roughness 
is  due  to  erosion,  while  the  large  scale  roughness  is  due  to  geophysical 
forces  that  result  in  hills  and  valleys,  or  as  in  the  case  of  the 
sea,  where  the  capillary  waves  are  dependent  on  surface  tension 
while  the  large  scale  rough  surface  is  generated  by  gravity  waves. 

For  the  general  case,  however,  one  cannot  assume  statistical 
independence  of  the  large  and  small  scale  surfaces. 

The  second  simplifying  assumption  that  was  made  was  that  the 

2  ' 

mean  square  slope  <?g  for  the  total  surface  was  approximately  equal 

2 

to  the  mean  square  slope  o^g  for  the  filtered  large  scale  surface. 

The  third  assumption  was  that  the  mean  square  height  of  the 
total  rough  surface  is  large  compared  to  a  wavelength,  and  the 


r 


surface  height  characteristic  function  for  the  total  surface  is 
negligibly  small  compared  to  unity. 

Finally,  the  physical  optics  approximation  for  the  cross 
polarized  backscatter  cross  section  is  zero  (Brown  1978).  As  a 
result,  the  cross  polarized  backscatter  cross  section  for  the 
filtered  surface  is  set  equal  to  zero  when  the  two- scale  model  is 
used.  However,  for  backscatter,  only  the  specular  points  on  the 
rough  surface  do  not  depolarize  the  incident  wave. 

In  Section  2.3  the  full  wave  approach  is  used  to  derive  a  unified 
formulation  for  the  like  and  cross  polarized  cross  sections  for  all 
angles  of  incidence.  These  solutions  are  compared  with  earlier 
solutions  based  on  a  two-scale  model  of  the  random  rough  surface 
(Bahar  et  al  1983).  Thus,  the  simplifying  assumptions,  that  are 
common  to  all  the  earlier  solutions  based  on  two- scale  models  of 
the  rough  surface,  are  carefully  examined.  It  is  shown  that  while 
the  full  wave  solutions  for  the  like  polarized  scattering  cross 
sections  based  on  the  two- scale  model  are  in  reasonably  good 
agreement  with  the  unified  full  wave  solutions,  the  two  solutions 
for  the  cross  polarized  cross  sections  differ  very  significantly. 

2.3  Application  of  the  Full  Wave  Solution  Without  Surface  Decomposition 

The  starting  point  for  this  analysis  is  the  full  wave  expression 


for  the  like  and  cross  polarized  scattering  cross  sections  of  the 
rough  surface  y  *»  h(x,z)  (Bahar  et  al  1983,  also  reported  in  Annual 
Technical  Report  March  1982-February  1983) 


k  __  D  P9 (n  ,n  |n)  , 

=  <SPQexp[iv  (h-h')]>  -  < - — - >  x(v  ) 

I-  ^  n»a  ^ 


in  which 


•exp[iv  x,  +  iv  z,]dx,dz, 
x  d  z  d  d  d 


r,  =  (x-x')a  +  (z-z')a  =  x ,a  +  za 

d  y  Z  dx  dz 


is  the  radius  vector  between  two  points  on  the  reference  plane  (x,z) 
The  vector  v  is 

v  =  k  (n^-n*)  *va+va+va  (2. 

o  xx  yy  zz 

where  kQ  is  the  free  space  wavenumber  for  the  electromagnetic  wave 
and  n*  and  n^  are  unit  vectors  in  the  directions  of  the  incident 
and  scattered  wave  normals  respectively.  An  exp(iut)  time  dependence 
is  assumed  in  this  work.  The  symbol  <  >  denotes  the  statistical 
average  and 

k^  k^  PQ  -  2 

““  <SP(*exp[iv  (h-h')]>  -  ~  R.~Snl  P2 (rif .n1  lns)p(hx»hz)dhxdhz 

•X2(vy,-vy)  =  IPQ(nf,ni)x2(vy,-vy)  (2. 

in  which 

n(h  ,h  )  is  the  unit  vector  normal  to  the  rough  surface 
X  z 

f(x,y,z)  -  y-h  (x,z)  *  0  (2. 


Vf  “  nlvfl  =  V(y-h(x,z))  =  (-h  a  +  a  -ha) 
11  x  x  y  z  z 

in  which  the  components  of  the  gradient  of  h(x,z) 


are  random  variables  and  p(h  ,h  )  is  the  probability  density  function 

X  z 

for  the  slopes  h^  and  h^.  The  expression  for  the  scattering  cross 
PQ 

sections  <o  > ,  (2.1)  is  valid  for  all  roughness  scales.  It  accounts 
for  shadowing  and 

P2<nf ,n*|n)  -  P2(n^,ni| ,ng)S(n^*n)S(-ni*n)+  (2.8) 

in  which  P2(n^,n*|n)  is  the  probability  that  a  point  on  the  rough 

surface  is  both  illuminated  and  visible  given  the  value  of  the 

slopes  at  the  point  (Smith  1967,  Sancer  1968)  and  P„(nf,ni|n  )  is  its 

z  s 

value  at  the  specular  points  where  the  unit  vector  n  is  given  by 

-f_-i 

n  -+•  n  «  —  -----  —  «  v/v  (2.9) 

s  |  — r  — i  | 

|n  -n  | 

The  arguments  of  the  unit  step  functions  S(-n^*n)  and  S(n^*n) 

vanish  at  points  of  the  rough  surface  where  the  incident  and 

scattered  waves  are  tangent  to  the  rough  surface.  Thus 

S(-n'I'*n  )  ®  1  and  S(n^*n  )  «  1  . 
s  s 

The  characteristic  and  joint  characteristic  functions  for  the 
surface  height  h  are  respectively  (Beckmann  and  Spizzichino  1963) 

x(vy)  -  <exp(ivyh)>  (2.10) 

and 

x2<vy»-vy)  *  <exp[ivy(h-h’)]>  .  (2.11) 

It  is  assumed  in  this  work  that  the  probability  density  function 
for  the  surface  height  is  jointly  Gaussian.  Thus 

2  2 

Sancer 's  expressions  for  P2  are  derived  for  kQ<h  >  >>1. 


(2.12) 


I X(Vy) 1 2  *  exp(-v2<h2>) 

and 

x2(vy.-vy)  =  exp(-v2<h2>  +  v2<hh’>)  (2.13) 

2 

where  <h  >  is  the  mean  square  height  and  <hh*>  is  the  surface 

''  PQ 

height  autocorrelation  function.  The  coefficients  D  x  depend  explic¬ 
itly  upon:  (1)  the  polarization  of  the  incident  wave  (second 
superscript  Q*V  -  vertical,  Q“H  -  horizontal)  and  (2)  the  polariza¬ 
tion  of  the  scattered  wave  (first  superscript;  P=V,H),  (3)  the  direction 
of  the  incident  and  scattered  wave  normals  n^  and  n^  respectively,  (4) 
the  unit  vector  n  normal  to  the  rough  surface  and  (5)  the  complex 
permeability  and  permittivity  of  the  medium  of  propagation  respec¬ 
tively  (Bahar  1981a).  On  deriving  (2.4),  it  is  assumed  that  the 
rough  surface  is  Gaussian  and  stationary,  thus  the  surface  height 
h  and  slopes  (hx,hz)  are  statistically  independent  (Brown  1978, 
Longuet-Higgins  1957).  It  is  also  assumed  that  for  distances 
| r^ |  less  than  the  surface  height  correlation  distance, 

£  ,n(h  ,h  )  =  n'Oi'jh').  It  has  been  shown  that  if  the  principal 
contributions  to  the  scattered  fields  come  from  specular  points 

on  the  rough  surface  (n  *  n  ) ,  (2,1)  reduces  to  the  physical  optics 

s 

solution  for  the  scattering  cross  section.  If,  however,  the 

roughness  scale  of  the  surface  is  small  compared  to  the  wavelength 
2  2 

(kQ<h  ><<1)  and  the  surface  slopes  h^  and  hz  are  very  small,  (2.1) 
reduces  to  the  perturbation  solution  for  the  scattering  cross 
sections  (Rice  1951).  Thus,  in  this  case  Bragg  scattering  is 


accounted  for  and  the  backscatter  cross  sections  for  grazing  angles 
are  strongly  dependent  on  polarization.  Recently  a  two-scale  model 
was  adopted  to  determine  the  corresponding  full  wave  solution  for 
the  scattering  cross  sections  (Bahar  and  Barrick  1983) .  To 
facilitate  the  application  of  the  two-scale  model  it  is  assumed 
that  the  small  scale  surface  hg  and  the  large  scale  filtered  surface 
hp  are  statistically  independent  (Valenzuela  1968;  Wright  1968; 

Brown  1978) .  This  assumption  is  reasonable  if  the  surfaces  hp  and 
hg  are  results  of  independent  processes  (Brown  1978)  as  for  example, 
when  the  small  scale  roughness  is  due  to  erosion  while  the  large 
scale  roughness  is  due  to  geophysical  forces  that  result  in  hills 
and  valleys  or  as  in  the  case  of  the  sea,  where  the  capillary  waves 
are  dependent  on  surface  tension  while  the  large  scale  surface  is 
generated  by  gravity  waves.  In  general,  however,  it  cannot  be 
assumed  that  the  large  and  small  scale  roughness  of  the  surface 
are  statistically  independent.  In  the  general  case,  if  the  two- 
scale  model  is  used  to  analyze  the  problem  it  would  be  necessary 
to  know  the  large  and  small  scale  surface  height  joint  probability 
density  function  for  two  adjacent  points  on  the  rough  surface  to 
determine  X£  (2.11)  alone. 

Since  the  full  wave  solutions  account  for  both  Bragg  scatter 
and  specular  point  scatter  in  a  unified,  self  consistent  manner, 
in  this  paper  solutions  for  (2.1)  are  developed  without  adopting 
a  two-scale  model  of  the  rough  surface. 
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It  has  been  noted  in  the  introduction  that  the  physical  optics 

approximation  for  the  cross  polarized  backscatter  cross  section 
PO 

is  zero  (<o  >  =  0  for  P^Q) (Brown  1978).  However,  even  the  large 

r 

scale  filtered  surface  will  depolarize  the  backscattered  field  at 

nonspecular  points  on  the  surface.  Therefore  the  present  analysis 

should  shed  more  light  on  the  evaluation  of  the  like  and  cross 

polarized  backscatter  cross  sections  and  the  suitability  of  the  two- 

scale  model  even  if  it  can  be  assumed  that  the  large  and  small 

scale  surfaces  are  statistically  independent. 

2  2  2 

Assuming  that  kQ<h  >  <<  1  and  |x|  <<  1,  the  scattering  cross 

section  (2.1)  can  be  expressed  as  follows: 


IP<^(nf  .n1)  fx2(vy,-vy) 


l  X  (Vy)  l  ] 


•exP[ivxxd  +  ivzzd]dxddzd 

s  I^Cn^n1)  Q  (n^,n*,R)  (2.14a) 

in  which  I  ^  is  defined  by  (2.4)  and 

oo 

Q  =  (X2  "  |x|2)exp[ivxxd  +  ivzzd]dxddzd  ,  (2.14b) 

— OO 

2 

is  the  two  dimensional  Fourier  transform  of  (x2  “  lx|  )  (2.12), 

(2.13).  It  therefore  depends  on  the  surface  height  correlation 
coefficient  R 

R  =  <hh' >/<h2>  . 


(2.15) 
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Using  the  notation  of  Rice  (1951),  the  surface  height  spectral 
density  function  W(vx>vz)  i s  related  to  the  two  dimensional  Fourier 
transform  of  the  surface  height  autocorrelation  function. 


W(v  ,v  ) 
x’  z 


_1 

7T2 


<hh’>  exp[ivxx^  +  ivzz^]dx^dzd 


(2.16a) 


and 


<hh'> 


(  W<vx’vz) 


exp[-iv  x,  -  iv  z,]dv  dv 
x  d  z  d  x  z 


(2.16b) 


Thus  assuming  that  the  rough  surface  is  Gaussian  and  stationary, 
to  compute  the  scattering  cross  sections  (2.14)  it  is  necessary  to 
prescribe  the  two  dimensional  slope  probability  density  function 


p(h  ,h  ),  (2.4)  and  the  surface  height  autocorrelation  function  or 
X  z 


its  Fourier  transform  (the  surface  height  spectral  density  function). 
Since  it  is  assumed  in  this  work  that  the  surface  is  isotropic, 

<hh'>  depends  only  on  the  distance  r^  *  |rd|  between  the  two  points 
(x,h,z)  and  (x'.h'.z')  on  the  rough  surface.  Thus 


<hh'>  =*  2ir 


W(vxz) 

4  ^o^Vxzrd^Vxz^Vxz 


(2.17) 


and 


<h  >  -  2tt 


W(v  ) 
xz 


v  dv 
xz  xz 


(2.18) 


in  which  J  is  the  Bessel  function  of  order  zero  and 
o 


xz 


2.2 
v  +  v 
X  z 


(2.19) 


Since 


J"(0)  =  [d2J  (v  r  )/dr2]  n  »  v2  /2 
o  o  xz  d  d  r  "0  xz 

d 


(2.20) 


it  follows  that 


R"(0)  *-o2/2<h2>  -  -o2/<h2>  »  -a2/<h2> 
S  x  z 


(2.21) 


Equation  (2.18)  for  the  mean  square  height  of  the  rough  surface 


k  +<  (k  -he  )  3(k  +k  ) 

c  c  c 


,2,  2  2.  3 

ok  (k  +k  ) 
c 

Thus  if  the  spectral  cut-off  point  k^  (Brown  1978)  is  much  larger 
than  k  (as  for  the  illustrative  example  (2,25)) 

<h2>  =  ~  =  87.9  cm2  (2.2 

6k 

The  surface  height  autocorrelation  function  <hh'>(2.17)  can  be 
expressed  in  closed  form  for  k^  ®.  Thus  the  surface  height 
correlation  coefficient  R(rd)  (2.15)  is  given  by  (Miller  et  al 
1972) 

R(rd)  =  [1  +  i  (xrd)2]  (Krd)K1(Krd)  -  (k^)2  K0(Krd)  (2.2 

in  which  Kq  and  are  the  modified  Bessel  functions  of  the  second 

kind  and  of  order  zero  and  one  respectively.  (See  Fig.  2.2) 

Since  k  >>  k  and  k  >  k  the  above  closed  form  expression  is  used 
c  c  o 

for  R  in  this  illustrative  example.  The  total  mean  square  slope 
of  the  rough  surface  is  obtained  on  substituting  (2.24)  into  (2.22) 


k  (6k  +15k  k  +  Ilk  ) 
c  _ c  _ c_ 

12(k2  +  k2)3 


(2.29). 


Thus  for  k  >>  < 
c 


T  .2  2  -i 

k  +  < 

B  yin  -£-5 - -  -  0.034 

K 


(2.31 


For  typical  sea  surfaces  the  relative  complex  dielectric  coefficient 

at  15  GHz  is  given  by  (Stogryn  1971) 

e  =  42  -  i39 
r 

The  slope  probability  density  function  p(h  ,h  )  is  assumed  to  be 

X  z 


Gaussian,  thus 


2  2 
.  n  +  n 

P<hx*hz)  "  ~  axp  ‘  ~~  2~ 
°S  °S 


(2.32 


In  Fig.  2.3  the  like  polarized  backscatter  cross  section  <o  >  is 

plotted  as  a  function  of  the  angle  of  incidence  0^  using  the 

expression  derived  in  Section  2.2  These  results  are  compared  with 

the  two-scale  full  wave  results  (Bahar  et  al  1983)  based  on  the 

choice  of  k^  (the  wavenumber  where  spectral  splitting  occurs) 

corresponding  to  8=1.  Both  results  yield  the  same  general 
W 

dependence  of  <o  >  on  the  angle  of  incidence.  The  small  difference 

in  level  is  primarily  due  to  the  fact  that  in  (2.3)  the  mean  square 
2 

slope  0^  of  the  total  (unfiltered)  surface  is  used,  (2.31),  while 

for  the  solution  based  on  the  two-scale  model  the  mean  square  slope 
2 

oFS  for  the  filtered  surface  hp  is  used  (Bahar  et  al  1983).  It 

should  be  noted  that  in  deriving  the  expressions  for  the  scattering 

cross  sections  based  on  the  two-scale  model,  it  was  assumed  that 
2  9 


2 


Thus  the  results  based  on  (2.3)  are  more  accurate. 


Backscatter  cross  sections  <ovv>  for  rough  surfaces  characterized  by  W(v  ,v  )  given  by 
(2.24)  with  B  =>  0.0046,  V  =  4.3  m/s  (a)  Two  scale  model  and  (A)  unifiedXfull  wave 
solution  (2,4).  Relative  complex  permittivity  is  e  »  42  -  i39,  A  *  2  cm. 


Furthermore,  on  deriving  the  solution  based  on  the  two-scale  model, 

the  quantity  X(v^)  (2.12)  is  assumed  to  be  negligible  compared 

to  X, (v  ,-v  )  (2.13)  for  r  <£  .  Since  4k^<h^>  *  3468  for  this 

illustrative  example,  the  resulting  approximation  is  very  good 

except  very  near  grazing  angles.  In  Fig.  2.4,  the  corresponding 

results  are  given  for  the  horizontally  polarized  backscatter 
HH 

cross  sections  <cr  >.  it  is  interesting  to  note  that  the  full  wave 

solution  (2.3)  yields  the  proper  polarization  dependence  of  the 

scattering  cross  sections  for  all  angles  of  incidence  without  use 

of  a  two-scale  model  since  it  accounts  for  specular  point  and  Bragg 

scattering  in  a  unified,  self-consistent  manner.  In  Fig.  2.5  the 

VH  HV 

cross  polarized  backscatter  cross  sections  <a  >  -  <a  >  are 

plotted  as  functions  of  the  angle  of  incidence.  Here  too,  both 

the  solutions  based  on  the  two-scale  model  as  well  as  the  solution 

derived  in  this  section  are  presented.  Unlike  the  solutions  for 

PP 

the  like  polarized  backscatter  cross  sections  <o  >  (P=V,H),  the 
solutions  for  the  cross  polarized  backscatter  cross  sections 
differ  significantly,  especially  near  normal  incidence  where  the 
difference  in  level  is  about  15db.  This  very  significant  difference 
is  due  to  the  fact  that  the  physical  optics  approximations  for  the 
cross  polarized  backscatter  cross  section  is  zero  (Brown  1978, 

Bahar  et  al  1983).  For  backscatter  the  surface  at  the  specular 
points  is  normal  to  the  incident  wave.  At  these  stationary  phase 
points  no  depolarization  occurs.  However,  since  depolarization 
occurs  at  the  nonspecular  points  of  the  filtered  surface,  the 
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physical  optics  approximations  for  the  cross  polarized  backscatter 

cross  section  is  not  valid.  It  is  interesting  to  note  that  for  the 

PP  PQ 

two-scale  model  at  normal  incidence  <a  >/<a  x>  =  47db  (P#Q) . 

PP  PO 

However,  using  the  full  wave  solution  (2.3),  <o  >/<a  y>  -  31db  (P^Q) . 

Consider  the  additional  illustrative  examples  at  X  band,  8.91  GHz 

A  =  3.367  cm,  k  =1.87  (cm)  ^ .  The  surface  height  spectral  density 
o  o 

function  W  is  given  by  (2.24).  The  wind  velocity  is  V  =  24  m/s  and 
the  relative  dielectric  coefficient  is  =  55  -  i37  (Stogryn  1971). 
Two  values  of  8  (2.24)  are  assumed 


B  ■  .0046 

<h2>  =  .853xl05  cm2 
o2  =  0.0498 


B  =  .0092 

<h2>  =  . 171x1 06  cm2 
o2  =  0.0997 


In  Figures  2.6,  2.7  and  2.8  the  backscatter  cross  sections  are 

evaluated  on  the  basis  of  the  unified  and  two-scale  models  for 

case  (a)  and  in  Figures  2.9,  2.10  and  2.11  the  cross  sections  are 

evaluated  for  case  (b) .  Daley  et  al  (1970)  published  extensive 

measurements  of  the  like  and  cross  polarized  backscatter  cross 

sections  at  microwave  frequencies  (Long  1975) .  From  his  data  at  X 

band  (8.91  GHz),  with  wind  speed  V  =  24  m/s, one  finds  that  for  an 

angle  of  incidence  6q  ■  15°,  <o^>/<o^>  =  20  db.  Using  the  unified 

PP  PQ 

full  wave  approach  <o  >/<o  >  =  21.5  db  for  case  (a)  and 

PP  PQ 

<o  >/<a  x>  -  18  db  for  case  (b) .  However,  using  the  two-scale  model 

PP  PO  PP  PQ 

<o  >/<a  *  36  db  for  case  (a)  and  <o  >/<o  x>  =  34  db  for  case  (b) . 


Fig.  2.8  Backscatter  cross  sections  <a  >  for  rough  surfaces  characterized  by  W(vx>vz) 

given  by  (2.24)  with  B  =  0.0046,  V  =  24.0  m/s  ( O)  Two  scale  model  and  (A) 
unified  full  wave  solution  (2.4).  Relative  complex  permittivity  is  e  =  55  -  i37, 
A  =  3.367  cm.  1 


2.9  Backscatter  cross  sections  <oV>  for  rough  surfaces  characterized  by  w(vx*vz)  given 
(2.24)  with  B  »  0.0092,  V  -  24.0  m/s  ( D)  Two  scale  model  and  (A)  unified  full  wave 
solution  (2.4).  Relative  complex  permittivity  is  e  =  55  -  i37,  X  =  3.367  cm. 


Thus  the  unified  full  wave  results  are  significantly  more  in  line 
with  published  experimental  results*  (Long  1975). 

2.5  Concluding  Remarks 

It  is  shown  in  this  section  that  two- scale  models  of  rough  surfaces 

can  be  adopted  to  obtain  solutions  for  the  like  polarized  backscatter 

cross  sections  that  are  in  reasonably  good  agreement  with  the  full  wave 

solutions  derived  in  this  section.  However,  the  two- scale  model  cannot 

be  used  to  evaluate  the  cross  polarized  backscatter  cross  sections.  The 

significant  differences  between  the  solutions  derived  in  this  section 

and  those  based  on  the  two-scale  models  are  primarily  due  to  the  fact 

that  the  physical  optics  approximation  for  the  cross  polarized  backscatter 

cross  section  (associated  with  the  large  scale  filtered  surface)  is  zero. 

For  backscatter,  the  specular  points  lie  on  portions  of  the  rough  surfaces 

that  are  perpendicular  to  the  incident  wave  normal  n^  (n  =  n^  =  -n*) . 

s 

At  these  specular  points,  the  backscattered  waves  are  not  depolarized. 
However,  at  nonspecular  points  of  the  rough  surface,  the  backscattered 
waves  are  depolarized  (Bahar  1981b).  Thus,  it  is  important  to  note  that 
even  if  a  surface  satisfies  the  radii  of  curvature  criteria  (associated 
with  the  Kirchhoff  approximations  for  the  surface  fields),  the  physical 
optics  approximations  for  the  scattered  fields  may  not  be  valid  unless 
for  the  given  incident  and  scatter  angles  specular  points  exist  on  the 
surface  and  significant  contributions  to  the  scattered  fields  come  from 
these  stationary  phase  (specular)  points  of  the  surface.  This  explains 

*See  also  Naval  Research  Laboratory  Report  7142,  Daley,  J.  C.,  W.  T. 

Davis  and  N,  R.  Mills,  "Radar  Sea  Return  in  High  Sea  States, 

September  25,  1970. 


why  the  physical  optics  approximations  for  the  like  polarized  backscattered 
cross  sections  are  not  suitabLe  for  grazing  angles  even  if  the  surface 
meets  the  radii  of  curvature  criteria  associated  with  the  Kirchhoff 
approximations . 

There  are  additional  important  reasons  for  preferring  to  use  the 
analysis  developed  in  this  section  over  those  that  are  based  on  two- 
scale  models  of  rough  surfaces.  Firstly,  if  the  two-scale  model  is 
used,  it  is  necessary  to  assume  that  the  large  and  small  scale  surfaces 
are  statistically  independent  (Brown  1978).  Secondly,  even  if  the 
assumption  of  statistical  independence  is  acceptable,  when  the  two- 
scale  model  is  used,  it  is  still  necessary  to  judiciously  specify  k^ 

(where  spectral  splitting  is  assumed  to  occur) .  These  problems  do  not 
arise  when  the  unified  full  wave  formulation  is  used  to  evaluate  the 
scattering  cross  sections. 
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